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Abstract 

Analytical treatment of time transfer problem for Earth-Satellite sys- 
tem is presented. The development was made in a complete relativistic 
framework. In accordance with modern clock precision and for low alti- 
tude orbits, we neglect the other bodies and consider only the \ Earth 
potential developed up to the J2 term in spherical harmonics. 



1 Introduction 

During the last few years, new technological developments increase significantly 
the time transfer accuracy using laser links. Further technological developments 
will go beyond the precision of the currently used relativistic theoretical descrip- 
tion generally used in data reduction programs. At present time, only the first 
order in the post newtonian metric is considered and the Earth potential is re- 
duced to the monopolar term, which corresponds to a spherical Earth (Blanchet 
& al. 2001). 

In foreseeable future laser links between ground based stations and low orbit 
satellites (T2L2, ACES) would reach precisions of order 10"^'* or even 10~^^ 
s. In this context, we show modelizing relativistic effects considering spherical 
Earth only is not sufhcient. Hence, we consider the other effects to select the 
ones which have to be taken into account at this level of accuracy. 

In sec. 2, we derive the orders of magnitude of the different relativistic 
effects, related to both Earth's potential and external bodies. We show, only 
the Earth's J2 potential term has to be included. In sec. 3, we give the analytical 
expressions of the trajectory and the corresponding propagation time. 



2 Selecting the relevant terms for low orbit satel- 
lites 

In general relativity, photons follow null geodesies of the space-time metric 5^^, 
which encodes the gravitational interaction. Following the IAU2000 convention 
(Soffel & al. 2003), we take the following form for the metric (one uses units 
such that G = c = 1) 
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ds^ = gy.^dxi^dx" = -(1 - 2U)dt^ + (1 + 2U)Sijdx'dx^ . (1) 



This metric is developed up to first order in the potential U, which represents 
the sum of the newtonian form of the potentials associated to each gravitational 
source. This form of the metric will turn out to be sufficient for the proposed 
application. The sources of gravitation are the Earth and external bodies. As 
usual, greek indices (a, /?,...) run from to 3 (space-time variables), while latin 
indices k,...) run from 1 to 3 (spatial variables). 

The Earth potential can be developed in spherical harmonics. The monopolar 
term (spherical Earth) writes 



where Me is Earth's mass. Beyond the monopolar term, the dominant term 
is the so called J2-term, which, in brief, corresponds to the Earth's oblateness 
contribution to the gravitational field. This term is of the order of 



where Re is Earth's radius. The other terms in the multipolar development 
are at least two orders of magnitude weaker. The (tidal) contribution of each 
external body (essentially the Moon and the Sun) to the potential is of the order 



where M and L are the mass and the distance of the involved body. 

UE,m has relative effects of the order of 10~® on the motion of a photon. 
Hence, for a ^ 1000 km, altitude satellite, the effect on the flying time of a 
photon linking a ground based station and this satellite can reach ~ 10~^^ or 
some 10~^^ s. Since J2 is of order 10~^, the contribution of this term is of order 
10""'^^ or some 10"^''' .s. The contribution of the other multipolar terms is then at 
best of order 10~^® s. It turns out that the contributions of external bodies are 
of the order of the J2 contribution for altitudes corresponding to geostationary 
satellites. Besides, for r < 17000 km, it would be meaningless taking external 
perturbations into account if harmonic terms beyond J2 are not included in the 
Earth's potential model. Second order terms in the metric (neglected in the 
metric presented) are of order 10~^^, hence induce time delay corrections of 
order 10~^° s. 

The precision which could be reached by time transfer experiments in a 
foreseeable future is of the order of lO"'^^ s. Hence, considering satellites at 
altitudes ~ 1000 km., it is necessary to include, beyond the Earth's monopolar 
term, the Earth's J2 contribution to the potential, but it is legitimate to neglect 
all the other contributions. 

3 Obtaining the time transfer 

As usual in rclativistic solar system experiments, we take the following form of 
the metric (constants G and c have been explicitly included): 



Ue,J2 ^ J2 



of 
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ds^ = g^.dx'^dx'' = _ (^1 _ ^ j c^dt" + 7^ j I dr |2 (2) 

which generahzes the ^ metric, in order to inchide viable alternative geometric 
gravity theories, like scalar-tensor theories (Will 1993). The deviation of those 
theories to general relativity is encoded in the 7 term which is unity in general 
relativity. Current experimental tests and astronomical observations show | 
7 — 1 I is less than some 10^^ (Will 2006). Hence, a possible deviation from 
general relativity should not have measurable effects in our problem. Since 
making computations in this more general context does not lead to dramatic 
complications, we present results in the general case, but it is always possible to 
set 7 = 1 if the general relativity case is considered. The light geodesic equation 
writes : 

{\h 1 

^ = '-k^k^d^g,. (3) 

with ka — gapk^, k" — A being an affine parameter along the geodesic. 

The wave vector fc" is isotropic which means k^ka = 0. At first order in the 
potential, the geodesic equation leads to: 



(^1) 
dfcp 

dA 



= 



(^1) 
dfc, 

dA 



(1 + 7) 



GMe 



-^Rs^di — 



— J2RE^di 
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(4) 
(5) 



HI) (0) (1/2) (1) 

where Q" = + Q" + is the development of Q" up to the first order 



in U of any quantity Q". Up to first order in U the solution reads 

(-1) 

a;° — Xq + f{x, n)X + c^^g{x, n) 

fix, n) = {l + 2^.hRl {A\ + B^}) 
g{x, n) = 

ln{-^^)+J2Rl[A n-x-,(^.-l^)+B n • aTo ( J, 



(6) 



)GMe 
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A, B, C are constant along the geodesic, Xq correspond to initial conditions and 
n is a normalized tri- vector (n ■ ft = 1) related to the components of the wave 

(^1) 

vector fc". We get from this formulae the function A (s;'^,^), from which we 

(-1) /(^i) ^ \ 
obtain ( A {x^,n) ) and find (writing, as usual, x'^ —t)- 



(^1) CAfj^ 
x' {t,n) = xl, + n'c{t-to) + ^nt-to,n) (7) 
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f (i-to,n) = (l+7) 
h\t~tQ,n) = ' 



r + X ■ n 
ro + xo ■ n 
1 1 



n ■ xq 



r 
9' 



n ■ xq 
1 

ro 



1 

'"0 



+ r(r-ro) 



with B\ A', 8*, F* constant along the geodesic. Obtaining the time trans- 

Hl) (-*!) (-.1) 

fer is then solving the equation a;* {t,n) —x^^^ (t), where x^^^ (t) = x^q + 

(-^1) (0) (1/2) (1) 

vlito)it~to) + ^{to){t-tof. We get T =T + T + T with T = t - to : 



(0) 



(1/2) 

T = 



^\f(x2o^~^hf (8) 

i/A (0) 
n •— T (9) 



(1) 1 
T= - 
2 



'(0) 



(0) 



(0) 

T 



(0) GAffi 
T T- 



(0) ^/(O) (0)^ 



(10) 



(0) 

ith = 



4 Conclusion 

As a consequence of the loss of spherical symmetry, a lot of new terms ap- 
pears from the J2 potential term ( see [7]). Possible next steps in order to 
upgrade our analysis could involve the next order in the metric development 
{goi — 0(3/2)) and/or the next terms in the spherical harmonic development 
( J3, C22, ■•■)• Thus, to go further, the next development could be based on a 3/2 
order metric with the earth potential developed up to the J2 term or more if 
needed. External potentials are not relevant for low altitude orbits and should 
be neglected in future studies. 

References 

[1] Blanchet, L., Salomon, C, Teyssandier, P., & Wolf, P. 2001, A&A, 370, 
320-329 

[2] Soffel, M., & al. 2003, Astronomical Journal, 126, 2687-2706 

[3] Will, C.1993, Theory and experiments in gravitational physics (Cambridge 
Univ. Press, Cambridge) 



[4] WiU, C. 2006, Living Rev. Relativity 



4 



